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Abstract 

In this paper, we consider a nj x multiple-input multiple-output (MIMO) channel 
subjected to block fading. Reliability (in terms of achieved diversity order) and rate 
(in number of symbols transmitted per channel use) are of interest in such channels. 
We propose a new precoding scheme which achieves both full diversity (ntn^th order 
diversity) as well as full rate {nt symbols per channel use) using partial channel state 
information at the transmitter (CSIT), applicable in MIMO systems including < nt 
asymmetric MIMO. The proposed scheme achieves full diversity and improved coding 
gain through an optimization over the choice of constellation sets. The optimization 
maximizes for our precoding scheme subject to an energy constraint. The scheme 
requires feedback of rat — 1 angle parameter values, compared to 2ntnr real coefficients 
in case of full CSIT. Error rate performance results for 3x1, 3x2, 4x1, 8x1 
precoded MIMO systems (with nt = 3, 3, 4, 8 symbols per channel use, respectively) 
show that the proposed precoding achieves 3rd, 6th, 4th and 8th order diversities, 
respectively. These performances are shown to be better than other precoding schemes 
in the literature; the better performance is due to the choice of the signal sets and the 
feedback angles in the proposed scheme. 



Keywovds — MIMO precoding, asymmetric MIMO, full rate, full diversity, partial CSIT, optimum con- 
stellation sets. 
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1 Introduction 



Multiple-input multiple-output (MIMO) techniques can achieve high data rates and spatial 
diversity in wireless communications over fading channels [1],[2]- Spatial multiplexing (V- 
BLAST) with rtt antennas at the transmitter achieves the full rate of rit symbols per channel 
use, but does not achieve transmit diversity. Space constraints in user terminals like mo- 
bile/portable receivers make asymmetric MIMO configuration with < Ut antennas at the 
receiver to be a preferred choice. Precoding techniques can improve performance through 
the use of channel state information at the transmitter (CSIT). Several precoding methods 
use CSIT and achieve diversity benefits, but compromise on the achieved rate; e.g., only 
one symbol per channel use is achieved in transmit beamforming [3]. It is desirable that 
precoding methods achieve high rates (preferably the full rate of rit symbols per channel use 
as in V-BLAST) and high diversity orders (preferably the full diversity order of ntnr)0 with 
partial CSIT. Our work reported in this paper addresses this problem. 

A vast body of research in the literature has addressed the problem of MIMO precoding. 
In particular, precoding using partial CSIT/limited feedback has been of interest because 
providing full CSIT (which refers to the full knowledge of all the channel gains between 
transmit and receive antennas) through feedback can be too expensive I1]-[Z]. 

Precoding in spatial multiplexing (V-BLAST) systems has been considered for achieving 
high rates and transmit diversity [9]-[T4j. Precoding schemes in [10]- [I2] incur some loss in 
rate. For example, the Grassmannian subspace packing based precoding in [10] does not allow 
simultaneous transmission of more than rif — 1 streams (i.e., achievable rate is < n^ — 1 symbols 
per channel use). The precoding scheme in [13] achieves full rate asymptoticalljl^ using D- 
BLAST architecture with partial CSIT, but achieves full diversity only for min(nf,nr) = 2 
and less diversity orders for other cases. Other recent works have proposed to improve the 
diversity gain of singular value decomposition (SVD) precodingfj by pairing good and bad 
subchannels, and jointly coding information across each pair of subchannels [15], [16]. In the 
schemes in [I5],[T6], though the achieved diversity orders are high, namely {nt — ^^ + l)(?T-r — 
^ + 1) where Ug < rij- is even number of streams, full diversity is not guaranteed. Also, 



^In this paper, 'full rate' is defined as rit symbols per channel use regardless of the number of receive 
antennas, and 'full diversity' is defined as ntUrth diversity. 

^The scheme in [T3] achieves a rate of 2T — 1 symbols per T channel uses in 2 x 2 system, which achieves 
close to full rate (of nt = 2) for large number of channel uses. 

■^SVD precoding achieves the full rate of rit symbols per channel use but its diversity order is only 1. 
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Prccoding 


Channel 
Knowledge 


Rate 


Diversity 


Grassmannian 
packing scheme [3] 


Partial CSIT 


Not full 
1 symbol/chl. use 


Full 

UtUr 


Grassmannian 
packing scheme [TP] 


Partial CSIT 


Not full 
< nt symbols/chl. use 


Not Full 


SVD prccoding 


Full CSIT 


Not full 
min(nt,nr) symbols/chl. use 


Not full 

1 


E-dmin Precoder in jl5j* 


Full CSIT 


Not full rate for rir < nt 

Hs symbols/chl. use 
Us < riri '^t ^ i^r, even 


Not full 
Us = imn(nt , nr) , Us even 


X-, Y-precoders in |1 61 * 


Full CSIT 


same as in [151 


same as in 1151 


Precoder in [T3]** 


Partial CSIT 


Asymptotically full 
2T — 1 symbols in T chl. uses 
for rit = rir = 2 


Not full for 
min(7if , nr) ^ 2 


Precoder in [14] 


Partial CSIT 
one angle 


Full 

nt symbols/chl. use 


Not full 


Proposed 


Partial CSIT 

nt — 1 angles 


Full 

nt symbols/chl. use 


Full 



Table 1: Comparison of required channel knowledge, achieved rates and diversity orders in different prc- 
coding schemes. 

* Fi-dmin and X-, Y-Codes in [TS] and [T^ differ mainly in complexity and applicability to higher-order QAM. 
** Precoding scheme in [T3] uses D-BLAST architecture which needs multiple channel uses to asymptotically 
achieve full rate. 



these schemes use full CSIT. In addition, the number of symbols transmitted per channel 
use Ug < Ur, and so for Ur < rit full rate is not achieved. The precoding scheme in [11] uses 
partial CSIT in the form of a single angle parameter feedback. Though this scheme achieves 
full rate and high diversity, it failed to achieve full diversity. A summary of the required 
channel knowledge, achieved rates and diversity orders in various precoding schemes are 
given in Table-dl 

In the above context, the significance of our contribution in this paper is that we propose 
a novel precoding scheme that achieves both full rate of rit symbols per channel use as well 
as full diversity of ntn^- using partial CSIT. The scheme requires the feedback of only — 1 
angle parameter values, compared to 2ntnr real coefficients in case of full CSIT. The proposed 
scheme achieves full diversity and coding gain through an optimization over the choice of 
constellation sets. The optimization maximizes for our precoding scheme subject to 
an energy constraint. We present codeword error performance results for nt = 3,4,8 and 
rij, = 1, 2. We analytically prove the full diversity {ntUrth diversity) of the proposed scheme, 
and simulation results are shown to validate the result. It is further shown that the proposed 
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scheme performs better than other existing schemes. This better performance is attributed 
to the choice of signal sets and feedback angles in the proposed scheme. 



Notation: Vectors are denoted by lower case boldface letters, and matrices are denoted by 
upper case boldface letters. (.)*, i-)^ , (.)^ and tr(.) denote conjugation, hermitian, transpose 
and trace operators, respectively. A^y will be used to denote the (x, y)th entry of matrix A. 

The rest of the paper is organized as follows. The system model is presented in Section 
[21 The proposed full rate, full diversity precoding scheme, and choice of constellation sets 
are presented in Section [3l Performance results and discussions are presented in Section HI 
Conclusions are given in Section [5l 



Consider a precoded V-BLAST system with nt antennas at the transmitter and rir antennas 
at the receiver. Let H G C"''^"'' denote the channel gain matrix, whose entries hpg are 
distributed i.i.d CJ\f{0, 1), = 1, ■ ■ ■ , n^, Vg = 1, ■ ■ ■ , n^. Let F denote the precoder matrix 
of size rit X rit, which is known to both the transmitter and receiver. For a given channel 
matrix H, the receiver computes rit — l number of real feedback parameters and sends them 
to the transmitter. Given the feedback parameters, the transmitter forms the precoding 
matrix F. Let x denote the rit x 1 complex data symbol vector of the form 



where Xj is transmitted by the ith antenna. Xj's will take values from a suitable complex 
constellation, the choice of which will be discussed later. The received signal vector, y G 
C"''^"^, at the receiver is given by 



2 System Model 



X = [ Xi X2 




(1) 



y = 



HFx + n. 



(2) 



where n G C' 



in,- X 1 



is the noise vector with its entries distributed as i.i 



.d. m{0,a^). 



3 Proposed Precoding Scheme 



Consider the precoded V-BLAST system described in Section [21 Let 




a a 



a ] , 



(3) 
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where the column vector a = [ai ■ ■ ■ druY ^ C"*^^ and |aj| = l,Vz = 1, ■ ■ ■ ,r2t. At 
high SNRs, symbol error probability decays as SNR""*"'', provided the received distance 
squared between each pair of codewords x^, {k I), denoted by dl is a Chi-square xlmnr 
distributed variable with 2ntnr degrees of freedom [T3]. In general, 

dl = llHFAxf, (4) 

where Ax = x^ — x/, A; 7^ /. Now, @ can be written as 

dli = tr(HFAxAx^F^H^) 

= tr(F^H^HFAxAx^). (5) 

It can be shown that the {p, q)th element of H^'^HF is given by 



nt rir 

fH^HF)„„ = (6) 



m=l 0=1 



It may be noted that all the elements of the pth row of the matrix H^HF are identical. 

Tit nt rtr 

^(F^H^HF),, = J2<{J2''-{T.^onhom)) 

n=l m=l 0=1 

nt nr nt nr 

= K\^{Y\^on\^) + Y<^rn[^J2Knhomy (7) 

n=m=l 0=1 nj^m o=l 

Moreover, it is clear that all the entries of the matrix F-'^H^HF are identical. So it amounts 
to 

(nt nr nt nr \ 

Yl l«nr(Xll^-l') + Yl <(^mY^onhon?)\B, (8) 

n=m=l 0=1 m,n,n>ni o=l / 

where Bpg = 1. In order to ensure that dli is distributed as a Chi-square xlmnr variable, we 
require that in ([8]), 

nt nr 

^( Y <Ci^n{Y^onhom)) = 0. (9) 

m,n,n>m o=l 

Let ai = \ai\e^'% YZ=i Knhom = I Eo=i /^L/iomle^"™, where j = v^. Then, with |ai| = 1, 
Vi = 1, ■ ■ ■ (jni) becomes 



nt nr 



Y I Y Kn^on] COs{9m - On + Onm) = 0. (10) 



m,n,n>m o=l 



For a given realization of H, the receiver needs to compute to satisfy (11 01) . Again, ffTOj) 
can be rewritten as 

nt n— 1 rir 



n=2 m=l 0=1 

Let the inner summation of (fTTj) with index m be equal to for each n. For n = 2, 



COs(6'i -02 + tt21; 



Let ^1 = 0. Hence, 62 = 021 — 7r/2. Next, for any n > 2, 

71—1 Tlr 



0. 



^71 I "-'nm 



m=l 0=1 

$Z I KriKm I (C0S(6'„ + COS 6^ + Sm{em + anm) slu On) 



n—l Tlr 



(12) 



(13) 



m=l o=l 



tan 



m=l 



n-1 



I E I Eo=l^on^om|sin( 
\ m=l 



I 



(14) 



Using (IT^ . 6'^ can be recursively calculated, given the values of 6i to 6'„_i. From ([S]), 
and ( fT4l) . after simplification, we have 



n=l o=l 



i=l 



where 



Ax = [ Axi Ax2 ■ ■ ■ Ax„,( ] . 



(15) 



(16) 



3.1 Choice of Constellation Sets 

From (lisp , it is clear that full diversity is guaranteed if, for each codeword pair x^, x^, /c 7^ / 
in the codebook, 



^Ax, ^ 0. 



(17) 



Let Xj, Vz = 1, ■ ■ ■ , 



TT-t, take values from some set Ci C ajZ[j] = {aiz : z G where 
^[j] = ^ + J^') ^ and w are integers, i.e., regular QAM constellations scaled by Oj. Clearly 
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then, Axi belongs to some other set Di C ajZ,[j] for each i. It may be noted that — Axj G -Dj. 
Now, let ai = 1 and choose 02 such that 

D,nD2 = {0}. (18) 

Again, define Di+ D2 = {Axi + Ax2 s.t VAxi G -Di, Ax2 G -D2}. Next, choose as such that 

D3n{Di + D2) = {0}. (19) 

Proceeding in this way, is chosen such that 

An(z^i + --- + A-i) = {0}. (20) 

This implies that for such choice of a^'s and Axi G Di, (ITTl) is satisfied. 

A trivial choice of a^'s could be y/pi, where p^'s are distinct prime numbers in Z. This shows 
that the solution set of Oj's satisfying ( I20l) is non-empty. For these values of a^'s, ( IT71) holds 
true for all codeword pairs in the codebook ensuring full diversity. Since the number of 
elements of Cj is finite, the cardinality of is also finite. This enables us to exhaustively 
enumerate all the elements of Di as functions of Oj. Next, we have to choose a^'s G C such 
that (I20I) is satisfied = 1, ■ ■ ■ Thus, there exist scaled QAM constellations for which 
the proposed precoding scheme guarantees full diversity. A parallel argument can be made 
for PSK constellation to ensure full diversity. 

Now, the d'^i^ parameter of a codebook determines its coding gain as given by the error 
probability expression at high SNR. Furthermore, d'^^^ is maximized if rnin{J2^^i ^^j) is 
maximized in ( fTSl) . Hence, after scaling the constellation sets Cj's by appropriate a^'s, we 
propose to rotate/scale each Ci by angles {(pi's) and scaling real numbers (fej's), ensuring that 
the average constellation power never exceeds the total transmit power constraint and also 
the condition (1201) holds true Vi = 1, ■ ■ ■ , n^. The corresponding min(^"^^ Axj) is computed 
over all Cj's. The optimum rotation angle {4>i,opt) and scaling real number ipi^opt) for each Ci 
are chosen by computer search for which min(S^^^Axj) is maximum. This fixes the choice of 
constellation sets Cj's that guarantee full transmit diversity as well maximizes subject 
to transmit power constraint. 
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3.2 A Simplified Approach to Choose Constellation Sets 



The treatment described in the previous subsection (Section 13. ip is feasible for constellations 
of small cardinality. The computational complexity involved in choosing large-sized constel- 
lation sets makes the choice cumbersome. In this subsection, we illustrate a suboptimal 
mechanism to choose large constellations. Our choice of F in (jS]) ensures that the effective 
transmitted vector Fx becomes 



It then implies that each of the nt antennas effectively transmits the symbol 'Y^^Li ^i- Now, 
for a given energy constraint, it is known that square QAM constellations (i.e., 4^^ sized 
constellation sets, where / is a positive integer) achieve better over PAM and PSK 
constellations. Hence, we propose to choose constellation sets for each Xi in such a way 
that "^^li Xi takes values from QAM constellation. Such a choice guarantees full diversity 
since for each codeword x transmitted the corresponding effective symbol Y17=i i^^ips to 
a unique point in some QAM constellation satisfying condition (fT7|) . 

3.3 Constellation Sets Obtained for rit = 3, 4, 8, 16 

Based on the treatment described in the above subsections (Sections 13. ![ 13. 2p . we carried 
out a computer search to obtain full diversity achieving constellation sets for various rtt. Let 
clQm = {o-q '■ q € QM,Ci G M}, where Qm denotes M-QAM constellation. The obtained 
constellation sets for nt = 3,4,8,16 are given in Table-2. Figures [D^a) and [D^b) show the 
plots of the obtained constellation sets for nj = 4 with of 1 bit /symbol and 2 bits/symbol, 
respectively. 

4 Results and Discussions 

We evaluated the codeword error rate (CER) performance of the proposed precoding scheme 
in asymmetric MIMO systems with < nt through simulations. In particular, we illustrate 
the CER performance of the proposed scheme for 3x1, 3x2, 4x1 and 8x1 full-rate MIMO 
as a function of the average received SNR per receive antenna. ML decoding is used. 

Plots of two performance measures, namely, i) CER and ii) pdf of the normalized ci^j^, are 




(21) 
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No. of Tx antennas 



nt = 3 



nt = 4 



Tit = 



nt = 16 



Modulation 
(1 bit /symbol) 



xi e {±1} 
X2 e {±j} 

X3 e {±0.675eJT} 



xi e {±1} 

2:^2 e {±j} 

X3 e {±i} 

a;4 e {±i?-} 



a;i G {±1} 
X2 e {±j} 

3^3 G {±5} 
0:4 G {±ij} 
3^5 G {±i} 
G {±ij} 
xre{±l} 
X8 G {±i_j} 



XI G {±1} 
X2 G {±j} 
X3 G {±i} 
0:4 G {±ij} 

3^5 G {±i} 

G {±ij} 
^7G{±i} 
G {±ij} 
e {± j^} 
xio G {±4j} 
xii G {±32} 

Xl2 G {±32^} 
X13 G {±i} 
^14 G {±64^} 
^15 e {±128} 
X16 G 



Modulation 
(2 bits/symbol) 



xi G Q4 
a;2 G iQ4 
a^a G jQi 



xi G Q4 

X2 G iQ4 
2^3 G |Q4 
2^4 G_g_Q4 



xi G Qi 

X2 G iQ4 
2:3 G |Q4 
2^4 G gQ4 

xe ^ 32 Q4 

2:7 G ^Q4 
Xs G 428^ 4 



G Q4 
a;2 G iQ4 
2:3 G |Q4 
a;4 G aQ4 
2:5 G YgQ4 

^6 ^ 52^4 

a;? G ^Q4 
Xs G 458Q4 

G 256 24 
XlO e 5pS4 

^11 € 10^4 24 

^12 S 20^24 

^13 € 40^24 

^14 G 81^2 Q4 

^15 € T6|84 24 



a;i6 G 



32768 



Q4 



Modulation 
(4 bits/symbol) 



a;i G Q16 

X2 G Y4Q16 
X3 G ^Ql6 



a;i G Q16 
a;2 G Y4Q16 

2^3 G 28 2l6 
X4 G 5fiQ l6 



2^1 G Q16 

a;2 G Y4Q16 
a;3 G ^gQie 
a;4 G ^Qi6 
G xp2i6 

2^6 G 224 2l6 
XT G 4i8Ql6 
^8 € 856 2l6 



a;i G Q16 

X2 G Y4Q16 

X3 G ^gQie 

•'E4 G ^Ql6 
2^5 G 4|2 2l6 
^6 G 224 2l6 
2:7 G 4|8Ql6 
^8 8^6 2l6 
2^9 G iT^Qie 
^10 £ 35p2l6 
^11 S 71p2l6 
^12 G T4p2l6 
^13 ^ 28!j)72 2l6 
^14 G 57|44 2l6 
■''l^ G 114^88 2l6 
^16 G 22fl37fi 2l6 



Table 2: Full diversity achieving constellation sets for nt = 3,4,8,16 with 1 bit/symbol, 2 
bits/symbol and 3 bits/symbol obtained by computer search. 



9 



(a) 1 bit/symbol, rit = 4 (b) 2 bits/symbol, rif = 4 

Figure 1: Plots of constellation sets for = 4 with 1 bit /symbol and 2 bits/symbol. 



shown in two separate sub-figures in each of Figs. [21 [3], |4]for 3x1,4x1 and 8x1 MIMO, 
respectively. For comparison purposes, we have plotted the CER performance and pdf 
of the full-rate precoding scheme in [T3]. We show the comparison only with [13] because 
both the proposed scheme and the scheme in [13] achieve the full rate of Ut symbols per 
channel use in < Ut asymmetric MIMO. Whereas, the other precoding schemes in Table-1 
lose rate and do not achieve full rate of rit symbols per channel use in rir < rit asymmetric 
MIMO. 

In Fig. [21 CER plots for the rates of 3 and 6 bits/channel use are shown. For 3 bits/cu, the 
constellation points used in the simulation for the proposed scheme are as given in the entries 
of rif = 3 and 1 bit/symbol modulation in Table-2. Likewise, for 6 bits/cu, the constellation 
points used are from the = 3 and 2 bits/symbol modulation in Table-2. Appropriate 
constellations from Table-2 are used for the plots of different rates shown in Figs. [3] and 
[H Since constellation optimization is not done in the precoding scheme in [Tl], we have 
used regular modulation alphabets (e.g., BPSK, 4-QAM) and matched the bits/cu in both 
schemes for fair comparison. For example, 6 bits/cu plot in Fig. [2] for the scheme in [H] is 
simulated using 4-QAM. 

From the CER plots of in Figs. [2] to [H we can see that the proposed precoding achieves 
better diversity orders compared to the precoding scheme in [Tl]. Indeed, as predicted by 
the analysis, the proposed scheme achieves the full diversity of ntn,. = 3,4,8 in Figs. [21 [3] 
and [H respectively. This can be verified by observing that pdfs of the ci^j^ in the proposed 
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scheme match with the theoretical xintn^ P*^^ (Chi square distribution with 2ntnr degrees of 
freedom). For example, in Fig. 121 the proposed scheme's pdf matches with that of the xl P'^f' 
verifying the achievability of ntrij. = 4th order diversity. Similar matches with Chi square 
distribution are observed in other figures as well. Thus, the analytical claim of full diversity 
of ritfir in the proposed precoding scheme is validated through simulations as well. It is 
further noted that the pdfs in the scheme in [13] do not match with those of the theoretical 
xlntur pdfs, indicating that the scheme in [M] does not achieve full diversity. Finally, Fig. 
[5] shows the performance of the proposed scheme in a 3 x 2 MIMO system with 3 and 6 
bits/cu, where the full ritrir = 6th order diversity is found to be achieved. 

5 Conclusion 

We presented a partial CSIT based precoding scheme, which achieved both full diversity 
(njn^th diversity) as well as full rate {rit symbols per channel use) with a feedback requirement 
of only fit — 1 real angular parameters, applicable in MIMO systems including asymmetric 
MIMO with Ur < rit. The full diversity was achieved by choosing optimized constellation 
sets. Through analysis and simulations we established the full diversity achievability in the 
proposed scheme. 
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Figure 2: Comparison of a) CER performance and b) pdf of the normalized ci^j„ of the 
proposed precoder with that in [2] for 3x1 MIMO system. 




Average Received SNR (dB) Normalized dmin 



Figure 3: Comparison of a) CER performance and b) pdf of the normahzed ci^j^ of the 
proposed precoder with that in [Tl] for 4x1 MIMO system. 
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Figure 4: Comparison of a) CER performance and b) pdf of the normalized of the 
proposed precoder with that in [13] for 8x1 MIMO system. 




Figure 5: CER performance and pdf of the normahzed of the proposed precoder in a 
3x2 MIMO system. 
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